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On  the  computation  of  space-time  correlations  by 
large-eddy  simulation 

By  Guo- Wei  Hef,  Meng  Wang  and  Sanjiva  K.  Lele 


1.  Motivation  and  objectives 

Sound  radiated  by  turbulent  flow  is  dependent  on  the  space-time  characteristics  of 
the  flow  field.  According  to  Lighthill’s  theory  (Lighthill  1952;  Proudman  1952),  acoustic 
power  spectra  in  the  far-fields  are  determined  by  the  two-time,  two-point  Eulerian  velocity 
correlations.  This  has  significant  implications  to  the  use  of  large  eddy  simulation  (LES) 
for  aerodynamic  noise  prediction.  The  existing  subgrid  scale  (SGS)  models  are  mostly 
constructed  to  predict  spatial  statistics  such  as  energy  spectra  (Meneveau  & Katz  2000). 
However,  it  is  not  clear  whether  these  models  can  lead  to  accurate  space-time  correlations, 
or  frequency  contents  at  individual  wavenumbers.  The  latter  information  is  essential  to 
accurate  noise  predictions,  because  for  a given  frequency,  only  the  spectral  element  of  the 
source  field  corresponding  to  the  acoustic  wavenumber  in  a given  direction  can  radiate  in 
that  direction  (Crighton  1975).  This  represents  a very  small  fraction  of  flow  energy,  and 
is  very  susceptible  to  numerical  and  modeling  errors.  Hence,  accurate  prediction  of  the 
space-time  correlation  is  a new  requirement  imposed  on  SGS  modeling  by  aeroacoustic 
applications. 

For  brevity,  we  henceforth  refer  to  the  two-time,  two-point  correlation  of  the  velocity 
field  simply  as  time  correlation.  It  describes  the  space-time  statistics  of  turbulent  flows. 
The  time  correlation  can  be  equivalently  expressed  by  a two-time  correlation  of  velocity 
Fourier  modes  in  spectral  space 


or  its  normalized  form 


C(k,  t ) = (iti(k,  t)uj(— k,  t + t)), 

(1.1) 

R(!.  _n  (ui(k,  t)ui(—k,  t + t)) 
(ui(k,t)ui(-k,t)} 

(1.2) 

Previous  research  (He,  Rubinstein  & Wang  2002)  compared  the  normalized  time  cor- 
relations or  correlation  coefficients  in  forced  isotropic  turbulence  calculated  by  direct 
numerical  simulation  (DNS)  and  LES  using  a standard  spectral  eddy  viscosity  model 
(Chollet  & Lesieur  1981).  The  comparisons  show  that  the  LES  overpredicts  decorrela- 
tion time  scales.  A following-up  study  (He,  Wang  & Lele  2002)  made  the  comparisons 
in  decaying  turbulence,  using  the  classic  Smagorinsky  model  (Smagorinsky  1963),  the 
dynamic  Smagorinsky  model  (Germano  et  al.  1991)  and  the  multiscale  LES  method 
(Hughes,  Mazzei  & Oberai  2001).  The  overpredictions  axe  still  observed.  These  numer- 
ical observations  motivated  the  present  analysis  to  investigate  the  overpredictions  and 
their  influence  on  sound  power  spectra.  The  analysis  will  be  carried  out  for  the  un- 
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normalized  time  correlations,  not  the  normalized  ones  as  before,  since  the  former  are  the 
ones  actually  used  in  the  computation  of  sound  power  spectra. 

The  analysis  starts  with  a generalization  of  Kraichnan’s  sweeping  hypothesis  (Kraich- 
nan  1964)  from  stationary  turbulence  to  decaying  turbulence.  This  involves  replacing 
a constant  convection  velocity  by  a time-dependent  one  in  a simple  kinematic  model. 
The  solution  of  the  kinematic  model  defines  a time-dependent  sweeping  velocity.  Kraich- 
nan’s sweeping  hypothesis  is  the  foundation  of  the  turbulence  theory  on  time  correlation. 
Kaneda  and  Gotoh  (1991)  and  Kaneda  (1993)  developed  the  Lagrangian  renormalization 
group  theory  and  the  Taylor  expansion  technique  for  time  correlations.  Rubinstein  and 
Zhou  (2000)  used  the  sweeping  hypothesis  to  formulate  the  scaling  law  of  sound  power 
spectra. 

The  present  analysis  on  time  correlations  will  shed  some  light  on  the  ability  of  LES 
to  predict  sound  power  spectra,  since  the  latter  can  be  analytically  expressed  as  an 
integration  of  time  correlations  using  the  quasi-normal  closure  assumption.  Previous 
evaluations  (e.  g.  Witkowska,  Juve  & Brasseur  1997;  Seror  et  al.  2001)  are  made  directly 
on  acoustic  fields.  These  calculations  unavoidably  have  to  cope  with  the  numerical  errors 
caused  by  the  truncation  of  the  source  region  (Crighton  1993;  Wang,  Lele  & Moin  1996). 
Instead,  we  will  discuss  the  influences  of  the  SGS  modeling  on  the  accuracy  of  sound 
prediction  through  an  analysis  of  time  correlations. 

In  this  paper,  we  consider  decaying  homogeneous  isotropic  turbulence.  Unlike  station- 
ary turbulence,  the  time  correlations  in  decaying  turbulence  are  dependent  on  both  time 
separations  and  starting  time.  Therefore,  two  different  starting  time  will  be  chosen,  one 
at  the  energy  propagation  stage  and  another  at  the  final  decaying  stage. 


2.  Main  results 

2.1.  Numerical  results 

A decaying  homogeneous  isotropic  turbulence  in  a cubic  box  of  side  27r  is  simulated 
by  DNS  with  grid  size  2563  and  LES  with  grid  size  643.  A standard  pseudo-spectral 
method  is  used,  in  which  spatial  differentiation  is  made  by  the  Fourier  spectral  method, 
time  advancement  is  made  by  a second-order  Adams-Bashforth  method  with  the  same 
time  steps  for  both  DNS  and  LES,  and  molecular  viscous  effects  are  accounted  for  by 
an  exponential  integrating  factor.  All  nonlinear  terms  are  dealiased  with  the  two-thirds 
rule. 

The  following  SGS  models  are  used  in  the  LES: 

(1)  The  spectral  eddy  viscosity  model:  we  use  the  Chollet-Lesieur  standard  form  for 
the  spectral  eddy  viscosity  (Chollet  & Lesieur  1981),  where  the  cutoff  energy  is  evaluated 
from  the  LES. 

(2)  The  Smagorinsky  model  (Smagorinsky  1963):  the  Smagorinsky  constant  is  Cs  = 
0.22  and  the  filter  width  is  set  equal  to  the  inverse  of  the  largest  effective  wave  number 
Ac  = 21. 

(3)  The  dynamic  Smagorinsky  model  (Germano  et  al.  1991):  the  Smagorinsky  coeffi- 
cients are  determined  by  the  Germano  identity.  The  grid  filter  width  is  k~x  and  the  test 
filter  width  is  taken  as  2fc~1. 

(4)  The  multi-scale  LES  method  (Hughes  et  al.  2001):  we  decompose  the  filtered  Navier- 
Stokes  equations  into  the  large  scale  equations  for  the  lower  one-half  Fourier  modes  and 
the  small  scale  equations  for  the  remaining  half  Fourier  modes.  The  dynamic  Smagorinsky 
model  is  only  applied  to  the  small  scale  equations. 
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model. 


The  initial  condition  for  DNS  is  an  isotropic  Gaussian  field  with  energy  spectrum 

E(k, 0)  oc  (fc/fc0)4exp[-2(fc/fc0)2],  (2.1) 

where  fco  = 4.68  is  the  wave  number  corresponding  to  the  peak  of  the  energy  spectrum. 
The  shape  of  the  energy  spectrum  excludes  the  effects  of  the  box  size.  The  initial  Reynolds 
number  based  on  Taylor’s  microscale  is  127.4.  The  initial  condition  for  LES  is  obtained 
by  filtering  the  initial  DNS  velocity  fields  with  filtering  wavenumber  kc  = 64/3  » 21. 
Therefore,  the  initial  LES  and  filtered  DNS  velocity  fields  are  exactly  the  same.  At  early 
stages,  the  LES  and  DNS  velocity  fields  are  highly  correlated  due  to  the  same  initial 
conditions.  Therefore,  the  time  correlations  of  the  LES  velocity  field  are  nearly  the  same 
as  those  of  the  DNS  field.  As  time  progresses,  the  LES  fields  become  decorrelated  from 
the  DNS  fields.  The  difference  in  time  correlations  between  the  LES  and  DNS  velocity 
fields  are  then  observed.  Therefore,  we  first  advanced  the  DNS  and  LES  velocities  in  time 
to  decorrelate  them  before  starting  to  calculate  the  time  correlations. 

The  energy  spectra  at  t = 0.5  and  t — 4.0  are  presented  in  Fig.  1.  Generally  speaking, 
the  LES  results  are  in  good  agreement  with  the  DNS  result  at  small  wavenumbers  but 
drop  faster  than  DNS  at  large  wavenumbers.  The  decay  of  the  total  resolved  energy  is 
presented  in  Fig.  2.  The  results  from  LES  with  all  SGS  models  follow  the  DNS  results 
with  some  deviations  throughout  the  entire  time  range.  They  exhibit  excessive  dissipation 
before  the  time  t — 1.5  (the  energy  propagation  range)  and  insufficient  dissipations  after 
t = 1.5  (the  final  decay  range).  Note  that  t = 0 does  not  correspond  to  the  start  time 
of  the  simulation.  Rather,  it  is  0.1  time  units  after  the  start,  and  hence  the  deviations 
among  DNS  and  LES  results  are  already  significant.  In  both  Figs.  1 and  2,  the  multiscale 
LES  results  are  the  best  of  all  models  and  the  classic  Smagorinsky  model  results  are  the 
worst.  The  dynamic  Smagorinsky  model  and  spectral  eddy  viscosity  model  yield  results 
in  the  middle  with  small  differences. 

Figure  3 plots  the  un-normalized  time  correlations  of  the  velocity  fields  from  the  DNS 
and  LES  for  wavenumbers  k = 5, 9, 13,  and  17,  spanning  a range  of  scales  from  the  integral 
scale  to  the  lower  end  of  the  resolved  scale.  The  starting  time  is  t — 0.5.  It  clearly 
shows  that  there  exist  discrepancies  between  the  LES  and  DNS  results,  and  that  the 
discrepancies  become  larger  with  increasing  wavenumbers.  Again,  the  classic  Smagorinsky 
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Figure  2.  Decay  of  total  resolved  energy. DNS; dynamic  Smagorinsky  model; 

multiscale  LES; Smagorinsky  model; spectral  eddy  viscosity  model. 
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Figure  3.  Time  correlation  C(k,r)  vs  time  lag  r with  starting  time  t = 0.5  for  (a)  k = 5, 

(b)  k — 9,  (c)  k = 13,  (d)  k — 17. DNS; dynamic  Smagorinsky  model; 

multiscale  LES; Smagorinsky  model; spectral  eddy  viscosity  model. 
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FIGURE  4.  Time  correlation  C(k,r ) vs  time  lag  r with  starting  time  t — 1.5  for  (a)  k = 5, 

(b)  k = 9,  (c)  k — 13,  (d)  k = 17. DNS; dynamic  Smagorinsky  model; 

multiscale  LES; Smagorinsky  model; spectral  eddy  viscosity  model. 

model  results  are  the  least  accurate  of  all  models,  the  multiscale  LES  method  is  the  most 
accurate,  and  the  spectral  eddy  viscosity  model  and  the  dynamic  Smagorinsky  model  lie 
in  between,  the  latter  being  slightly  better  than  the  former. 

Figure  4 plots  the  same  time  correlations  as  in  Fig.  3 but  with  a different  starting  time 
t = 1.5.  The  discrepancies  observed  are  qualitatively  the  same  as  in  the  t = 0.5  case, 
which  indicates  that  the  SGS  modeling  errors  equally  affect  the  time  correlations  in  the 
final  decay  range. 

In  summary,  it  is  observed  in  decaying  isotropic  turbulence  that  discrepancies  exist 
between  the  un-normalized  time  correlations  calculated  from  DNS  and  those  from  the 
LES.  These  discrepancies  will  be  analysed  in  the  next  section. 

2.2.  Analysis  of  the  numerical  results 

The  analysis  is  based  on  the  generalized  sweeping  hypothesis  for  decaying  turbulence. 
In  the  sweeping  hypothesis  for  stationary  isotropic  turbulence,  the  convection  velocity 
is  constant  (Kraichnan  1964).  However,  in  decaying  turbulence,  the  convection  velocity 
varies  with  time.  A generalization  can  be  made  by  introducing  a time-dependent  convec- 
tion velocity,  which  evolves  slowly  relative  to  the  time  scales  of  velocity  fluctuations. 
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Consider  a fluctuating  velocity  Fourier  mode  u(k,  t)  convected  by  a bulk  velocity  v(t). 
We  assume  that  the  wavenumbers  k of  the  fluctuating  velocity  are  sufficiently  large.  The 
flow  scales  associated  with  these  wavenumbers  are  small,  over  which  the  convection  ve- 
locity is  spatially  uniform  and  relatively  large  in  magnitudes.  In  this  case,  the  convection 
effect  is  dominant.  The  governing  equation  for  the  fluctuating  velocity  modes  is  therefore 


<9u(k,  t ) 
dt 


+ i[k  • v(t)]u(k,  f)  -0, 


which  yields 

u(k,  t + r)  = u(k,  t)  exp  i k • v(s)ds  ) . 
Then,  the  time  correlation  can  be  expressed  by 


(2.2) 


(2.3) 


<u(k,t  + r)u(-k,f)) 

/ I rt+r  pt+T  ' \ 

= (u(k,  t)u(— k,  t))  exp  ( — -k1 2  J J (v(s  )v(s  ))ds  ds  ).  (2.4) 

In  the  derivation  of  (2.4),  the  convection  velocity  v(t)  is  assumed  to  be  Gaussian  and 
independent  of  the  velocity  u(x,  t)  at  the  starting  time  t.  These  assumptions  can  be 
justified  by  the  near-Gaussianity  of  the  large-scale  velocity  and  its  initial  independence 
of  the  small-scale  velocity.  By  introducing  a sweeping  velocity 

1 rt+r  ft+T 

y2(t>T)  = ^2  J J (v(s')v(s  ))ds  ds" , (2.5) 

we  obtain  a general  expression  of  time  correlation  similar  to  the  one  in  stationary  tur- 
bulence 

(u(k,  t + t )u(— k,  t))  = (u(k,  t)u(— k,  t))  exp  ^-^fe2F2(t,r)r2^  . (2.6) 

The  calculation  of  the  sweeping  velocity  (2.5)  can  be  further  simplified  by  assuming 
the  following  form  of  the  bulk  velocity  correlation 

(v(s)v(s"))  = <v2(s'))exp  (— A|s'  - s"|)  , (2.7) 

where  A-1  is  a correlation  time  scale.  Substituting  (2.7)  into  (2.5),  we  find 

V2(t,r)  = -j-  J (v2(s  ))A_1  ^2  - exp[-A(s,  - t)]  - exp[-A(t  + r - ds  . (2.8) 

In  isotropic  turbulence,  the  bulk  velocity  is  determined  by  large  scale  motions.  Hence, 
its  decorrelation  time  scale  A-1  is  much  larger  than  those  of  velocity  fluctuation  modes 
considered  here.  Since  the  time  separation  t of  interest  is  within  the  decorrelation  time 
scales  of  the  velocity  fluctuations,  we  have  At  1.  Using  Taylor  series  expansion  with 
respect  to  At  and  ignoring  the  second  and  higher  order  terms  in  (2.8),  we  obtain 

1 /’t+T 

V2(t,r)  = — / (v2(s  ))ds' . 

T Jt 


(2.9) 
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Figure  5.  Normalized  time  correlation  R(k,r)  vs  (a)  un-normalized  and  (b)  normalized  time 

lag  for  different  modes  computed  using  DNS.  k = 5; K = 9; k = 13; 

k = 17;  □ k = 30;  a fc  = 40;  v k = 50;  > k = 60;  < k = 70;  o k = 80. 


Note  that  the  bulk  velocity  is  associated  with  the  energy-containing  motions,  and  its 
variance  (v2(t))  is  the  total  energy  of  the  instantaneous  velocity  field.  Since  the  energy 
decay  is  relatively  small  over  the  decorrelation  time  scale,  the  sweeping  velocity  can  be 
easily  calculated  by  simple  approximations  V2(t,  r)  = (v2(t  + r))  or  V2(t,  r)  = [(v2(f))  + 
(v2(t  + r))]/2.  The  former  is  used  in  the  following  discussions. 

Figure  5 plots  the  normalized  time  correlations  R{k,r)  from  DNS  for  wavenumbers 
k = 5,9, 13, 17, 30, 40, 50, 60, 70  and  80,  where  the  correlations  are  normalized  by  the 
instantaneous  energy  spectra.  The  time  separation  is  un-normalized  in  Fig.  5(a)  and 
normalized  by  the  scale-dependent  similarity  variable  Vk  in  Fig.  5(b).  The  latter  figure 
clearly  exhibits  that,  with  time  normalization,  all  curves  for  different  modes  collapse. 
This  verifies  the  validity  of  the  generalized  sweeping  velocity  in  decaying  turbulence. 

Equation  (2.6)  indicates  that  the  normalized  time  correlations  are  solely  determined 
by  the  sweeping  velocities.  In  the  present  LES,  the  sweeping  velocities  are  smaller  than 
that  of  DNS  values  because  of  the  reduced  total  energy.  Therefore,  the  time  correlations 
in  LES  decay  more  slowly  than  the  ones  in  DNS.  That  is  to  say,  the  LES  overpredicts 
the  decorrelation  time  scales  compared  to  DNS.  Figure  6 plots  the  normalized  time 
correlations  from  the  DNS  and  LES  with  respect  to  the  un-normalized  time  for  the  modes 
k = 5, 9, 13, 17.  It  confirms  that  the  time  correlations  from  LES  decay  more  slowly  than 
the  ones  from  DNS.  Again,  the  multiscale  LES  method  is  the  most  accurate  and  the 
classic  Smagorinsky  model  is  the  least  accurate  of  the  all  models  tested.  The  dynamic 
Smagorinsky  model  and  the  spectral  eddy  viscosity  model  are  in  the  middle. 

Equation  (2.6)  also  indicates  that  if  the  time  separation  is  normalized  by  Vk,  the  un- 
normalized  time  correlations  are  solely  determined  by  the  instantaneous  energy  spectra. 
Figure  7 plots  the  un-normalized  correlations  vs  the  normalized  time  separation.  It  shows 
that  the  LES  underestimates  the  magnitudes  of  time  correlations  relative  to  the  DNS 
results.  The  underestimation  obviously  becomes  more  significant  as  the  wavenumber 
increases,  which  is  consistent  with  the  more  severe  drops  of  the  LES  energy  spectra  at 
high  wavenumbers.  Again,  the  relative  performance  of  the  SGS  models  in  terms  of  the 
magnitudes  of  time  correlations  is  the  same  as  before. 

In  conclusion,  the  discrepancies  between  the  time  correlations  computed  using  DNS 
and  LES  consist  of  two  parts:  the  decorrelation  time  scale  and  magnitude.  The  errors  in 
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Figure  6.  Normalized  time  correlation  R(k,r ) vs  time  lag  r with  starting  time  t = 0.5  for 

(a)  k = 5,  (b)  k = 9,  (c)  k = 13,  (d)  k - 17. DNS; dynamic  Smagorinsky  model; 

multiscale  LES; Smagorinsky  model; spectral  eddy  viscosity  model. 


decorrelation  scales  are  induced  by  the  instantaneous  sweeping  velocity,  and  the  errors 
in  magnitudes  are  induced  by  the  instantaneous  energy  spectra.  In  relative  terms,  the 
errors  in  decorrelation  time  scales  are  less  significant  than  those  in  magnitudes.  However, 
they  must  not  be  ignored  since  the  sound  power  spectra  are  sensitive  to  the  decorrelation 
time  scale,  (see  discussions  in  the  next  section).  Note  that  the  sweeping  velocity  used  in 
our  analysis  is  the  r.m.s.  of  velocity  fluctuations,  or  the  square  root  of  the  total  energy. 
Thus,  an  accurate  prediction  of  the  instantaneous  energy  spectra  is  most  critical  to  the 
accurate  computation  of  the  time  correlations. 

2.3.  Discussion 

It  is  difficult  to  determine  the  influence  of  SGS  modeling  on  acoustic  power  spectra 
through  numerical  evaluation  in  general,  since  it  involves  the  truncation  of  the  source 
region  (Wang  et  al.  1996).  To  avoid  this  problem,  we  use  an  analytical  expression  of  acous- 
tic power  spectra  based  on  Lighthill’s  theory  and  the  quasi-normal  closure  assumption. 
The  analytical  expression  is  only  valid  for  stationary  turbulence.  However,  reasonable 
inferences  can  be  drawn  for  decaying  turbulence  through  this  analysis. 


Space-time  correlations  by  large-eddy  simulation 


327 


Figure  7.  Un-normalized  time  correlation  C(fc,r)  vs  the  normalized  time  lag  r with  starting 

time  t = 0.5  for  (a)  k = 5,  (b)  fc  = 9,  (c)  k = 13,  (d)  k --  17. DNS; dynamic 

Smagorinsky  model; multiscale  LES; Smagorinsky  model; spectral  eddy 

viscosity  model. 


According  to  Lighthill’s  theory  (Lighthill  1952),  the  acoustic  pressure  in  a far-field 
position  x is  given  by 


p(IC’,)=4 (2'10) 

where  Ty(y,t)  = pui(y,t)uj(y,t)  is  the  Lighthill  stress  tensor,  $7  the  source  region,  p 
the  mean  far-field  density,  c the  speed  of  sound  in  the  far-field,  and  y a position  vector 
in  the  source  field.  Based  on  this  equation,  Kraichnan  calculated  the  far-field  acoustic 
power  spectral  density  (Kraichnan  1953) 


w 


p M = ( \niniTa 


(2.11) 


where  n is  the  unit  vector  in  the  direction  of  the  point  x,  u the  frequency  of  the  radiated 
sound,  and  (k,  w)  the  space-time  Fourier  transformation  of  the  Lighthill  stress  tensor. 
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(a) V = 1.0; V = 0.95; V = 0.90.  (b) V = 1.0  with  full  spectra; 

V = 0.95  and  kc  = 25; V = 0.95  and  kc  = 13. 


Using  the  quasi-normal  hypothesis,  the  acoustic  power  spectral  density  function  can  be 
written  in  the  form  (Rubinstein  & Zhou  2002) 

P(“)  = f',7^/+  /L  ^!(*.Oexp(-»r)dr  (2.12) 

In  the  following  discussion,  the  normalized  time  correlation  R(k,r)  is  taken  as  the 
exponential  form 


R{k,  t)  = exp  (~fcW^  . (2.13) 

and  the  energy  spectrum  E(k)  assumes  the  von  Karman  form 

E(k)  = Ce2/3*0-5/3(fc/*o)4[l  + (k/k0)2}-17'6k-2,  (2.14) 

where  fco  = 5 defines  the  peak  of  energy  spectrum.  The  exponential  form  and  the  von 
Karman  spectrum  are  the  appropriate  approximations  to  the  time  correlations  and  the 
energy  spectra,  respectively,  in  our  numerical  simulations. 

With  the  substitution  of  (2.13)  into  (2.12),  the  non-dimensionalized  sound  power  spec- 
tra are  given  by 


PtM  = ~-pM^  j*°° (-4^2)  (215) 

where  M = Kj/c  is  the  Mach  number  and  Vo  = u>o/&o-  fco  is  the  inverse  integral  length 
scale  and  ojq  the  inverse  integral  time  scale. 

The  influences  of  the  overpredicted  decorrelation  scales  on  acoustic  power  spectra  can 
be  seen  in  Figure  8(a),  where  the  sound  power  spectra  are  evaluated  according  to  (2.15) 
with  the  sweeping  velocities  V equal  to  1.0,  0.95  and  0.9.  The  small  variations,  up  to 
10%,  of  the  sweeping  velocities  cause  significant  reductions  of  the  sound  power  spectra 
at  higher  frequencies.  This  illustrates  the  sensitivity  of  the  acoustic  power  spectra  to  the 
sweeping  velocities. 
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The  sweeping  velocity  induced  errors  can  be  compounded  by  the  truncation  of  the 
energy  spectra  at  high  wavenumbers,  corresponding  to  unresolved  scales  in  LES.  To  test 
this  effect,  the  energy  spectrum  is  truncated  ( E(k ) set  to  zero)  for  either  k > 25  or 
k > 13.  Figure  8(b)  plots  the  acoustic  power  spectra  for  the  sweeping  velocity  V = 0.95 
with  the  above  truncated  energy  spectra.  It  shows  that  the  acoustic  power  spectra  drop 
considerably  at  moderate  to  high  frequencies,  and  the  spectral  peaks  are  shifted  towards 
left  to  lower  frequencies. 


3.  Conclusions  and  future  work 

Numerical  comparisons  in  decaying  isotropic  turbulence  suggest  that  there  exist  dis- 
crepancies in  time  correlations  evaluated  by  DNS  and  LES  using  eddy-viscosity-type 
SGS  models.  This  is  consistent  with  the  previous  observations  in  forced  isotropic  turbu- 
lence. Therefore,  forcing  is  not  the  main  cause  of  the  discrepancies.  Comparisons  among 
different  SGS  models  in  the  LES  also  indicate  that  the  model  choice  affects  the  time  cor- 
relations in  the  LES.  The  multi-scale  LES  method  using  the  dynamic  Smagorinsky  model 
on  the  small  scale  equation  is  the  most  accurate  of  the  all  models,  the  classic  Smagorin- 
sky model  is  the  least  accurate  and  the  dynamic  Smagorinsky  model  and  spectral  eddy 
viscosity  model  give  intermediate  results  with  small  differences. 

The  generalized  sweeping  hypothesis  implies  that  time  correlations  in  decaying  isotropic 
turbulence  are  mainly  determined  by  the  instantaneous  energy  spectra  and  sweeping  ve- 
locities. The  analysis  based  on  the  sweeping  hypothesis  explains  the  discrepancies  in  our 
numerical  simulations:  the  LES  overpredicts  the  decorrelation  time  scales  because  the 
sweeping  velocities  are  smaller  than  the  DNS  values,  and  underpredicts  the  magnitudes 
of  time  correlations  because  the  energy  spectrum  levels  are  lower  than  the  DNS  ones. 
Since  the  sweeping  velocity  is  determined  by  the  energy  spectra,  one  concludes  that  an 
accurate  prediction  of  the  instantaneous  energy  spectra  guarantees  the  accuracy  of  time 
correlations. 

An  analytical  expression  of  sound  power  spectra  based  on  LighthilPs  theory  and  the 
quasi-normal  closure  assumption  suggests  that  the  sound  power  spectra  are  sensitive  to 
errors  in  time  correlations.  Small  errors  in  time  correlations  can  cause  significant  errors 
in  the  sound  power  spectra,  which  exhibit  a sizable  drop  at  moderate  to  high  frequencies 
accompanied  by  a shift  of  the  peaks  to  lower  frequencies. 

Based  on  the  above  analysis,  two  possible  ways  to  improve  the  acoustic  power  spectrum 
predictions  can  be  considered.  The  first  is  to  construct  better  SGS  models  to  improve 
the  LES  accuracy  for  time  correlations.  The  second  is  to  remedy  the  temporal  statistics 
of  the  Lighthill  stress  tensor  in  order  to  “recover”  the  contribution  from  the  unresolved 
scales  in  LES  to  time  correlations. 
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